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Abstract. We consider a generalized Lieb-Liniger model, describing a one- 
dimensional Bose gas with all its conservation laws appearing in the density matrix. 
This will be the case for the generalized Gibbs ensemble, or when the conserved charges 
are added to the Hamiltonian. The finite-size corrections are calculated for the energy 
spectrum. Large-distance asymptotics of correlation functions are then determined 
using methods from conformal field theory. 

PACS numbers: 02.30.1k, 05.30.Jp 
1. Introduction 

There has recently been much interest in so called generalized Gibbs ensembles (GGEs), 
where the density matrix of a system is given by 

Pgge = Zch E e-E.**-, (1) 

generalizing the usual Gibbs ensemble for generic (non-integrable) systems where only 
the Hamiltonian and particle number operator appear in the exponent. Instead, for an 
integrable system, {Q n } is the complete set of local conserved charges, with {(3 n } the 
generalized inverse temperatures (Lagrange multipliers) and Zqge — Tr e~^™ /3n< ^' 1 the 
generalized partition function pp. For such quantum ensembles, expectation values of 
observables are obtained as 

{6) GGE = Tr pggeO, (2) 

whereas the time-evolution is governed by the usual Schrodinger Hamiltonian. 

Recent experimental advances with quantum non-equilibrium dynamics of cold 
atoms P], where the GGE was proposed [3] to describe local large-time behavior, has 
spurred a great interest in possible equilibrium ensembles for integrable systems (see e.g. 
Refs. jU |5j [6] and references therein, and in particular Refs. [3 EJ [9J [10] for treatments of 
the one- dimensional Bose gas). But effects from higher conservation laws in integrable 
systems have also been studied longer ago in the context of competing interactions in 
spin chains (TTJ [12] . 
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It is therefore a timely question to now ask in general what sort of effects one 
can anticipate when incorporating higher conservation laws into the density matrix for 
an integrable model. In this paper we study finite-size effects and obtain conformal 
dimensions for the one-dimensional Bose gas (Lieb-Liniger model), when governed by a 
density matrix of the form (TjQ). In particular, using the exact Bethe Ansatz solution we 
calculate the finite-size corrections for the energy and momentum. Comparing this to 
the expressions given by conformal field theory, we are able to obtain the large-distance 
asymptotics of correlation functions. 



2. Lieb-Liniger model 

The Lieb-Liniger model describes a one-dimensional Bose gas with point-like interaction 
through the Hamiltonian 



H 



dx 



d x ^(x)d x ^f(x) + c^ t (a;)^ t (a;)^(x)^( 



x) 



(3) 



with coupling constant c > and Bose fields \1/ with equal-time commutation relations 
^(x), &(y) = 5(x - y) and [V(x), V(y)) = &(x), &(y) = . The model is solved 
through the Bethe Ansatz [13], [HI [15], yielding eigenstates |{Aj}) given in terms of 
rapidities Xj satisfying the Bethe equations 



N 



-n Ai 



A fc + ic 



J 



N, 



(4) 



for a system of N particles in a box of length L with periodic boundary conditions. The 
Bethe equations (J4j) can equivalently be written as 



N 



27m 7 -, 



J 



1 N, 



(5) 



k=l 



with rij integer when iV odd and half-integer when iV even, and 9(X) 
i In [(ic + A)/(zc — A)]. For non-zero wave function, rij ^ nk when j ^ k. 
The conserved charges Q n have eigenvalues Q n given by 



Qr 



N 



(6) 



the three lowest being particle number = Qq, momentum P = Qi, and energy 
E = Q 2 . Explicit expressions for some of the higher conserved charges Q n can be found 
in Ref. [TB] . Now let us consider a generalized Hamiltonian Hq where all the conserved 
charges Q n have been added to the Lieb-Liniger Hamiltonian ([3]), 

H G = H + J2 bnQn (7) 
with coefficients b n . Extremizing the entropy gives a density matrix of the form (TjQ) 



Pg 



v-i p -m G 

Zj g e . 



(8) 



Finite-size effects from higher conservation laws for the one- dimensional Bose gas 3 

with (3 = 1/T the inverse temperature and Zq the generalized partition function. Hence 
we can analyze both the situation of a system described by a generalized Gibbs ensemble 
(JI|) as well as a system where the Schrodinger Hamiltonian itself is given by Hq in Eq. (JTj). 
The eigenvalues E({Xj}) of the generalized Hamiltonian Hq are given by 

N 

= £ £ o(A,) (9) 
where the bare one-particle energy £"o(A) is given by the polynomial function 

oo 

e (A) = £ b n X n (10) 

n=0 

with b 2 — 1, and &o — — ^ the chemical potential. It was shown in Ref. [17] that the 
generalized Yang- Yang thermodynamic Bethe Ansatz (TBA) equation becomes 

£ ( A ) + ^ / d f iK(\,n)Hl + e- pe V>) = e (\), (11) 

Z7Tp J-oo 

and that a solution exists provided that Eq(\) is bounded from below and 
limA->±oo Eq(X) = +oo. The kernel K{X,p) is given by 

K(X, p) = 9'(X -fi) = 2c/ (c 2 + (A - y) 2 ), (12) 

and the equilibrium particle distribution function p(X) by 

p(A) = i?(A) ( J- + J- r dp if (A, /i)p(/i)) (13) 

VZ7T Z7T J-oo / 

with the Fermi weight 

"' A > = rrW- (14) 

This can also be written as p(A) = i9(X)pt(X), where 

ft (A) - ^ J™JiM${y)K{\p) Pt {y) = ± (15) 
is the density of states. Similarly, the so-called dressed charge Z(X) is defined by 

Z(X) - ±- r dp$(p)K(X,p)Z(p) = 1, (16) 

Z7T J-oo 

reflecting that for the Lieb-Liniger model Z(X) = 2irp t (X). 

3. Finite-size corrections at zero temperature 

3.1. Energy 

Let us now investigate the finite-size corrections to the generalized energy (|9~1), focusing 
on the zero-temperature limit. Then the set of numbers rij in Eq. (jSJ) that minimizes 
the generalized energy is such that all Xj with e(Xj) < are occupied and the rest 
empty of particles. In the thermodynamic limit, this gives Fermi points qf wherever 
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e(qf) = 0, with a filled Fermi sea i for < A < qf . The ground state energy in the 
thermodynamic limit can then be written as 

/oo rqi~ 
dXp(X)e (X) =LY, r dXp t (X)e (X), (17) 
-oo ■ J a. 



with finite-size corrections coming from the replacement of the sum in Eq. ([9]) with an 
integral. These are obtained using the Euler-Maclaurin formula in complete analogy 
with the usual case [T8| [19], yielding for the leading corrections 

Eo = LY, f! d\ Pt (\)s (\)-^-Y,Ml (18) 

i % i,v 

where i is the Fermi sea index and v = ± the right/left index. Here vf is the Fermi 
velocity at the Fermi point q i , given by 

1 de 



vt 



(19) 



Now we investigate the finite-size corrections of low-energy excited states, using 
the techniques found in Refs. [201 Ell 1221 [23J [2U [25]. First, we expand the energy 
E({qf}) = Eq + 8E to second order around the ground-state energy Eq in Eq. ( TT81) . 



with ^g? 1 the change in Fermi momentum g^ with respect to the ground state. Now, 

S = ±M(?f )eo(gf ) + ^E <*A %^e„(A) = ±M(gf ),(gf ) (21) 

where e is the dressed energy ( TTTT) . Since the ground state minimizes E, one has 

£ (g±) = 0. (22) 
Then Eq. ( f20|) becomes 

^ = § E v&fM-jL [PtiliXQi)} (23) 
which is given by 

& E = § E^(<?rw) (^n 2 = tt^eki [p t (gr)^n 2 (24) 

Now we will express Jgf 1 in terms of the quantum numbers AiVj and A.Dj, defined 
as the changes in 

Ni = Lj_ dpp t (p), (25) 

a=^( r - r)dpp t (p), (26) 
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compared to the ground state. Then 



3 J J ■> 



The Jacobian is found from 



9Ni 
dqf 



±L P t{qt 



Lpt{qt 



8ij ± 



dpg{p\qf 



\j_ 3 oo ~j+) d v9(p\qt) 



dqf 

where the function g(\\qf) is defined through 

An k q k An 

Eq. ( l30i) is equivalent to g{\\qf) 
to consider the vectors 

( Nt \ 



(27) 

(28) 
(29) 

(30) 



dpt(X)/dq i /pt{qf)- For the general case we need 



n 



D 1 

V ; J 



V 



fqi\ 

qt 

q% 

V = / 



(31) 



M l3 



1 



dni 



so that we can write 

Ani = L^2M ij p t (p j )5p j , ^ Lpt{pj)dp . 

with drill dp j given by Eqs. ff28]) - f|29]) . and hence obtain p t (q i )5qf from 

1 



Pt(pi)$Pi 



(32) 



(33) 



Inserting this into Eq. (I2lj) . and including the numbers N { and of particle- hole 
excitations at q~ and qf respectively, gives the general expression 

2NT + (E( M-1 ki A ^) ■ (34) 



5E 



-E 



This expression is simplified when the model has parity symmetry, e.g. around 
A = so that only even powers of A appears in Eq. fflOl) and hence eo(A) = e (— A). If this 
is the case, we can consider the new pairs of Fermi points qf = ±qj (qi > q 2 > ... > 0), 
with a sea j of either particles or holes between them, and Fermi velocities = Vi. 
Now define 

m 



Da 



dpp t (p), 



9j 



-1j 



1j J 



dpp t (p), 



(35) 
(36) 
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i.e. ANj being the number of particles/holes added to Fermi sea j of particles/holes, 
and dj = ADj/2 the number jumping from qj to —(jj. Let us now write 



dqf 
dDi 



±L Pt ( gj )2 



-Hi 



= ±Lp t { gj )2Z, 



(37) 
(38) 



with g{\\qf) given by Eq. fl30|) . By expressing the matrix elements of Z as Z, 
o~ij + f^dpgijiff), with g^ff) = (-l)^ 1 ^/^-), one finds % = where f(A) is 

the dressed charge matrix defined through 



1 f°° 

tijW - T7~E / df* K*k(\ n)€kj(n) = Sij, 

Z7T ; J-oo 



(39) 



so that here the kernel matrix element is i£y(A, /i) = K(\,fi) when — qj < // < qj, 
otherwise zero. Now from 

= - E I ** dX IdkiW - gu(-X)} z kj , (40) 

A= 9l k J -°° 



Z%j $ij 



it follows that 



Y 



8ij + 



dfigij(fi) = 5ij - [Zki - hi] (Z = (Z 



(41) 



Hence Y = (Z T y 1 . Putting this into Eqs. (J2ZD an d (J24J), finally gives 



a: 



(42) 



where 



\ fc k I 



(43) 



Obviously, when the bare single-particle dispersion in Eq. (FIU|) is given by £q(A) 



A , there is just a single Fermi sea, between A = — q and A = q, and Eqs. ( 1421) and (T43 
reduce to those for the usual Lieb-Liniger model, 



5P = ^ 
L 



A+ + A" 



A d 



iV + 



1 / AN 



2 \2Z{q) 



±Z(q)d 



(44) 



3. ,2. Momentum 

Since the momentum P = Y,i \ — x Si as is easily seen from Eq. (J5J), it is trivially 
obtained as 



P = F„ + £ AJV, -k i d i + '^- {N+ -Nr + d.AN,} 



L 



(45) 



for the excited states, where Po is the ground-state momentum. In the thermodynamic 
limit p = (g+ + q^)/2 and ki = qf — q~ . In case of parity symmetry Pj = and 
kj = 2 qj . 
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4. Conformal dimensions and correlation functions 



It is now clear that the finite-size corrections in Eqs. (1181) . (1421) and (1451) can be written 
on the form 

7T 



E - E (L -> oo) = 



A+ + Al 



h.o.c. 



p-p« = Z [p, AJVi - 



h.o.c. 



(46) 
(47) 
(48) 



where Tt.o.c. denotes higher-order corrections. This tells us that the low-energy physics is 
given by a sum of conformal field theories [26j EZ] , where A^ are the scaling dimensions 
of the scaling operators of the theories and Cj = 1 the central charges. We will now use 
this to obtain the large-distance asymptotics of the equal-time correlation functions [21J. 
The field correlation function 



($(x)&(0)) = Tr p G ^(x)^(0), 

for which the total number of added particles is AiV = ^-(A/V^-i — AN 2 j] 
then have the asymptotic form 



<*(x)*t(o))~£A(Q*) e 



*E,^'M' X -2E,( A J + + A 7) 



(49) 
1, should 

(50) 



where Qq, = {{ANj}, {dj}, {Nj } \ AN = 1} is the set of quantum numbers for the 
excitations, and A(Q^) the amplitudes. The leading term is given by 

(*(x)* f (0)) ~ x~ a , (51) 

with a the smallest sum of scaling dimensions 2J2j{A^ + Aj) given the constraint 

AiV = 1, i.e. a is the smallest of the numbers (j2j(Z~ 1 )jkJ /2. Similarly, density- 
density correlators are obtained with AiV = 0, 

(j(x)j(o)) - (j(o)) 2 ~ £ MO*) 



x -2V.(A++A-) 



(52) 



where j(x) = ^(x)^(x) and Q d = {{AiV,}, {dj}, {Nf} \ AN = 0}. The leading terms 
are 



(j(x)j(O)) - (j(0)> 2 ~ A x x 2 + A 2 cos{d k k k x)x 



(53) 



± 



where the first term comes from the processes with d k = and one single number N t 
equal to one, and the second term is from the two processes with the d k = ±1 giving 
smallest 9, i.e. with 6 the smallest of the numbers 2 Z k ^j ■ 

Importantly, the finite-temperature correlation functions for mixed states are 



obtained by the standard conformal mapping z(w) 



,2-irTw/vi 



z = x — ivit, so that in the formulas above 



x 



-2A d 



-»■ {(irT/v j )/amh[irTx/v j ]} 



2A d 



in the complex plane 



(54) 
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5. Discussion 

We have obtained the finite-size corrections for the energy and momentum of a 
one-dimensional Bose gas with delta-function interaction (Lieb-Liniger model) when 
additional conservation laws are present in the density matrix, or the conserved charges 
are added to the Hamiltonian. The results show that the low-energy physics is described 
by a sum of conformal field theories each with central charge c = 1, where each may 
have its own specific speed of sound. The picture is most clear when the bare dispersion 
is parity symmetric. In this case we derived the asymptotic behavior of the correlation 
functions using standard arguments. 

The finite temperature mapping ([Ml) provides a possible connection to present 
studies of quantum non-equilibrium dynamics. In the standard setup, an isolated system 
in a pure state is time-evolved by, but not an eigenstate of, the usual Schrddinger 
Hamiltonian. The density matrix of a subsystem will however be in a mixed state, 
presumably approaching the form (TjQ) with a finite effective temperature. It remains to 
be seen whether the new types of correlation effects in the generalized model studied 
here may appear in such systems. 

It is interesting to note that the generalized one-dimensional Bose gas with a 
dispersion relation with many Fermi points for the bare particles gives similar types 
of finite-size effects as in well-studied examples of multicomponent (a.k.a. nested) Bethe 
Ansatz solvable models [29] . such as integrable quantum spin chains [TH [121 I2S EQl 
EH E21 E31 EI] an d the one-dimensional Hubbard model [2U EE] , even though the one- 
dimensional Bose gas only contains a single species of particles. 
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